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Resumen 


Se verifica la relación, básicamente lineal, que existe entre la longitud de una 
melga y el gasto óptimo de riego, que es el gasto que se debe aplicar para obtener 
un valor máximo en el coeficiente de uniformidad, que produzca valores elevados 
de las eficiencias de aplicación y de requerimiento de riego. Dicha verificación 
se realiza mediante la aplicación de un modelo hidrodinámico completo, 
empleando métodos numéricos, que hace uso de las ecuaciones de Saint-Venant 
para la descripción del flujo del agua sobre el suelo, acopladas internamente con 
la ecuación de Richards, que permite modelar el flujo del agua en el suelo. De 
acuerdo con los resultados obtenidos mediante la aplicación del modelo para diez 
tipos de suelo diferentes, se incluyen cuadros de diseño operativo del riego por 


melgas. 
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Introducción 


El objetivo del diseño de riego por melgas es 
conseguir la aplicación de la lámina de riego 
requerida por el cultivo de la manera más 
uniforme posible, para alcanzar una eficiencia 
de aplicación alta. 

El diseño del riego consiste en determinar 
el gasto Óptimo y el tiempo durante el cual se 
aplica dicho gasto en la cabecera de la melga 
para lograr la mayor uniformidad de aplicación 
posible, para una longitud y ancho de melga 
específicos, y un tipo de suelo en particular. 

Rendón et al. (1997) han verificado que el 
gasto óptimo de diseño guarda una proporción 
lineal con la longitud de la melga en la cual 
debe aplicarse. El resultado es obtenido 
haciendo uso de un modelo formado por las 


ecuaciones de Lewis y Milne (1938), para 
describir el flujo del agua sobre la superficie del 
suelo, y de Green y Ampt (1911), para describir 
el flujo del agua en el suelo. Sin embargo, es 
necesario tener en cuenta que las formas de las 
ecuaciones utilizadas por los autores citados 
son de las más simplificadas en el contexto de 
la modelación del riego por gravedad, motivo 
por el cual en este trabajo se verifica la relación 
de proporcionalidad entre el gasto óptimo y 
la longitud de la melga, haciendo uso de las 
ecuaciones de Saint-Venant y Richards. 

En este trabajo se utiliza un modelo hidro- 
dinámico completo desarrollado por los autores, 
el cual hace uso de las ecuaciones de Saint- 
Venant para la descripción del flujo del agua 
sobre el suelo, acopladas internamente con la 
ecuación de Richards, que permite modelar 
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el flujo del agua en el suelo, con el objetivo 
de verificar la relación de proporcionalidad 
que existe entre la longitud de una melga y el 
gasto óptimo de riego, que es el gasto que se 
debe aplicar para obtener un valor máximo en 
el coeficiente de uniformidad, que produzca 
valores elevados de las eficiencias de aplicación 
y de requerimiento de riego. De acuerdo con 
los resultados obtenidos mediante la aplicación 
del modelo para diez tipos de suelo diferentes, 
se presentan cuadros de diseño operativo del 
riego por melgas. 


Flujo del agua sobre la superficie del 
suelo 


El flujo del agua con una superficie libre 
se modela con las ecuaciones de Saint- 
Venant que resultan de la aplicación de las 
leyes de conservación de masa y cantidad 
de movimiento. En una melga, la relación 
entre su ancho y el tirante de agua permite 
considerar las ecuaciones correspondientes al 
escurrimiento sobre una superficie de ancho 
infinito (Woolhiser, 1975). La ecuación de 
continuidad se escribe como: 


óh, 9,0 _ ad) 
ot ox ot 


La ecuación de momentum se escribe en la 
forma recomendada por Saucedo et al. (2011): 


2 
E ls 1 esos, 


hot mox 1% y Jox 

q ol 
+p2"=0 2 
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donde q(x,t) = Ulx, t)h(x,t) es el gasto por 
unidad de ancho de melga [L?T*]; x, la coor- 
denada espacial en la dirección principal del 
movimiento del agua en la melga [Ll]; t, el 
tiempo [T]; U, la velocidad media; h, el tirante 
de agua [L]; J, la pendiente topográfica de la 


melga [LL"]; J, la pendiente de fricción [LL”); 
V, = 01 / 0t, el flujo de infiltración [LT”], es decir, 
el volumen de agua infiltrado en la unidad de 
tiempo por unidad de ancho y por unidad de 
longitud dela melga; I, la lámina infiltrada [L]; g, 
la aceleración gravitacional [LT?]; el parámetro 
adimensional f = 1 — a, siendo a =1 - U,,/U, 
donde UI, es la proyección en la dirección del 
movimiento de la velocidad de salida de la 
masa de agua debido a la infiltración. 

La relación entre las variables hidráulicas q 
y h con la pendiente de fricción, denominada 
ley de resistencia hidráulica, es adoptada de 
acuerdo con Fuentes et al. (2004), haciendo uso 
de una ley potencial de resistencia: 


¿ont 
q= E) (3) 
v 


donde v es el coeficiente de viscosidad cine- 
mática del agua [L?T”] y k es un factor adi- 
mensional. 

Para establecer las condiciones inicial y de 
frontera que deben sujetar a las ecuaciones de 
Saint-Venant para modelar el riego por melgas, 
se tiene en cuenta una melga cerrada, para 
evitar la pérdida de masa fuera del dominio de 
riego. Para la fase de avance se tiene: 


q(x,0)=0 y  h(x,0)=0 (4) 


q(0,1)=q.,  a(xp.t)=0 y h(x,,t)=0 (5) 
donde x/(t) es la posición del frente de onda 
para el tiempo t y q, es el gasto de aporte 
en la cabecera de la melga. Para la fase de 
almacenamiento: 

g(0,t)=49,,  q(L,t)=0 (6) 
donde L es la longitud de la melga y q, el gasto 
de aporte en la cabecera de la melga. Para la 
fase de consumo: 


g(0,t)=0, g(L,t)=0 (7) 


Saucedo et al., Gasto óptimo en riego por melgas 


Para cerrar el sistema es necesario conocer 
la forma en que evoluciona en el tiempo la 
lámina infiltrada en toda posición sobre la 
melga, es decir, la ley de infiltración, la cual 
es proporcionada con base en la solución de la 
ecuación de Richards (1931). 


Flujo del agua en el suelo 


La descripción del proceso de infiltración del 
agua en el suelo puede efectuarse teniendo 
como base la ecuación de Richards (1931), 
que resulta de la combinación de la ecuación 
de continuidad con el campo de velocidades 
calculado conforme a la ley de Darcy y que en 
su forma tridimensional —sin tener en cuenta 
la extracción de agua por las plantas— se 
escribe como: 


ly. ON 
EV [KV] , (8) 


donde y es el potencial de presión del agua 
en el suelo, expresado como la altura de una 
columna equivalente de agua [L] (positivo en 
la zona saturada y negativo en la zona no satu- 
rada del suelo); C(yw) = d8/dy es denominada 
la capacidad específica de humedad del suelo; 
9 = O(1p) es el volumen de agua por unidad de 
volumen de suelo o contenido volumétrico de 
agua [L*L*], y es una función de y conocida 
como curva característica de humedad o 
curva de retención del agua; K = K(y) es la 
conductividad hidráulica [LT"], que en un 
suelo parcialmente saturado es una función del 
potencial de presión; el potencial gravitacio- 
nal es asimilado a la coordenada espacial 
z, Orientada positivamente hacia abajo [L); 
V = (0/0x, 0/dy, 0/02) es el operador gradiente; 
x, y son las otras dos coordenadas espaciales 
[L], y t es el tiempo [T]. 

El riego es un fenómeno que se efectúa en 
tres dimensiones, lo cual hace necesario que el 
flujo del agua en el suelo sea descrito utilizando 
la ecuación (8), sin embargo, es conveniente 
aceptar la hipótesis de que el fenómeno se 
efectúa en planos paralelos al desarrollo de 


la melga; entonces es posible utilizar la forma 
bidimensional de la ecuación de Richards: 


o) 
OX 


ay yl o 
Ama ROA 


Ko 221) (9) 
OZ 


Habida cuenta de los resultados de Saucedo 
et al. (2006) sobre el grado de aproximación 
de la hipótesis del tiempo de contacto en el 
riego por melgas, es posible simplificar más 
la aproximación, empleando la ecuación de 
Richards en su forma unidimensional: 


wy 0 y 
C(yw)—==—|K(yw)| —=- (10) 
me 2 E ] 
que debe resolverse sobre el dominio de 
solución, definido por una columna de suelo. 
Como condición inicial para la solución de 
la ecuación de Richards unidimensional, se 
debe especificar la distribución de las presiones 
en el espacio: 
y=Y.(2) (11) 
En la superficie del suelo se impone una 
condición de frontera tipo Dirichlet, con el 
potencial de presión igual al tirante de agua 
calculado mediante las ecuaciones de Saint- 
Venant: 
y=h, z=0, t>s0 (12) 
En la frontera inferior se tiene en cuenta una 
condición de gradiente unitario: 


Ay) _ 3 
0Z 


ñ z=P, t>0 (13) 


Donde P es la profundidad de la columna 
de suelo, que debe ser mayor que la máxima 
posición que el frente de humedecimiento 
alcanza durante el tiempo en que transcurre el 
riego. 
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La solución de la ecuación de Richards 
hace indispensable representar las propieda- 
des hidrodinámicas del suelo, expresando el 
potencial de presión (yw) como una función 
del contenido volumétrico de agua (0) y la 
conductividad hidráulica K como una función 
de 6. Como es señalado por Fuentes et al. 
(1992), la combinación de las características 
hidrodinámicas de Fujita (1952) y Parlange et 
al. (1982) es conveniente en estudios teóricos, 
como la construcción de soluciones analíticas 
exactas; en estudios experimentales puede ser 
más conveniente la combinación de la curva de 
retención propuesta por van Genuchten (1980), 
considerando la restricción de Burdine (1953), 
con la curva de conductividad hidráulica 
propuesta por Brooks y Corey (1964), debido a 
que satisfacen las propiedades integrales de la 
infiltración y a la facilidad para identificar sus 
parámetros. 

La curva de retención propuesta por 
van Genuchten (1980) se escribe como: 


[e(1) -e, ]/[8, 8, ]= Le(w/9)"] > donde V, 


es un valor característico de la presión del 


agua en el suelo; m y n son dos parámetros de 
forma empíricos relacionados por la restricción 
de Burdine (1953): m = 1 - 2/n, con 0<m< 
1 y n > 2; 6, es el contenido volumétrico de 
agua a saturación efectiva del suelo, y 0, es 
el contenido volumétrico de agua residual. 
La conductividad hidráulica propuesta por 
Brooks y Corey (1964) se representa como 
K(0)=K, [(e-0,)/(8, -0,)]”, donde n es un 
parámetro de forma empírico y positivo. 

La lámina infiltrada se calcula para cada 
posición x, donde se resuelven numéricamente 
las ecuaciones de Saint-Venant: 


1(x,,t)= 


o— 


[0(x,,2,t) -0(x,,2,0)]az (14) 


donde F es la posición del frente de hume- 
decimiento y 0(x,z,0) es el contenido volu- 
métrico de agua inicial. 


Solución numérica 


Solución numérica de las ecuaciones de Saint- 
Venant 


Para el cálculo del gasto óptimo, se hace uso de 
un esquema numérico para las cuatro fases del 
riego por melgas desarrollado por Saucedo et 
al. (2011). La disposición de las celdas de cálculo 
para las fases de avance, almacenamiento, 
consumo y recesión, se muestra en la figura 1. 
La forma discreta de la ecuación de continuidad 
para la fase de avance se escribe como: 


(03, +(1-0)q ]t-(x, -x,)[0(1,+1,) 
+(1-0)X%, +1,)] 

-[009, +(1-09)9 y ]98+(x, —X; )[0(1, + 1,) 
+(1-0)H(%,, +1,,)] (15) 
—[ 4h, +(1- 9)h, +91, +(1- 9)1, |(x, -x,) 
+[01,+(1- 9), +91, +(1-0)1,,](%,, =x,)=0 


para las fases de almacenamiento, consumo y 
recesión se escribe como: 


[3, +10), Jo1-[09, +10), 
—[9,+(1-9)M, +91, +(1-0)1, [x, -x,) (16) 


+[01; +(1- 4), +01; +(1- 9), -x,)=0 


La ecuación de momentum guarda la misma 
forma discreta para las cuatro fases del riego: 


23h(g, -q,)ot+ [gh -7?)(a, —h Jl 

+H(x, -x,)[09,+(0 0,09, - (1-0), | 
+90 x,-x,) [0], +(1-0)],-J,] (17) 
+potqh(x, -x,)ol, +(1-0)1, -oL,, 


-(1-0)1,|=0 


En las ecuaciones (15) a (17), Ót es el paso de 
tiempo, w y q son factores de peso en espacio 
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Figura 1. Disposición de las celdas de cálculo para las fases 


de avance, almacenamiento, consumo y recesión. 


y tiempo, y las demás variables se definen de 
acuerdo con la figura 1. 

El cálculo de los coeficientes se realiza con 
base en los valores pertenecientes al nivel de 
tiempo anterior q = (1 -— p)%,, + 99, h=1- 4) 
h, + dh, (Saucedo et al., 2011). Se introducen 
las variables pequeñas (Strelkoff y Katopodes, 
1977): 6h, 9q,, 8h, y 0q,, tales que, para las celdas 
interiores, se tiene: h, =h, + 0h, h,=h,, + 0h, q,= 
q,+0q, y q,=3,, + 0q,. Se sustituyen en las formas 
discretas de las ecuaciones de continuidad 


m 


y momentum, y se obtiene un sistema de 
ecuaciones algebraicas, cuya solución permite 
avanzar en el tiempo los valores de las variables 
hidráulicas. 


Solución numérica de la ecuación de Richards 


La ecuación de Richards unidimensional (10) 
es discretizada en el espacio, utilizando ele- 


mentos finitos y, en el tiempo, mediante un 
esquema implícito en diferencias finitas. El 
procedimiento se encuentra bien documentado 
en la literatura y puede consultarse, por 
ejemplo, en Saucedo et al. (2002). 


El gasto óptimo 


Para el diseño del riego por melgas se pueden 
utilizar dos métodos (Rendón et al., 1997): 
pruebas de riego y modelos físico-matemáticos. 
Las pruebas de riego consisten en hacer 
experimentos de campo, en donde las variables 
experimentales son el gasto de riego, la longitud 
de la melga y el tipo de suelo. Con base en los 
resultados de los experimentos, para cada tipo 
de suelo se obtienen combinaciones de gasto 
de riego unitario y longitud de la melga, que 
hacen posible aplicar con un alto grado de 
uniformidad una lámina de riego específica. 
Las pruebas de riego tienen la desventaja 
de requerir un tiempo considerable para su 
ejecución, son costosas y debe realizarse una 
cantidad elevada de pruebas para determinar 
el gasto de riego para el que se obtiene la mayor 
uniformidad en el riego; por estas razones, son 
recomendadas sólo para calibrar los métodos 
de diseño. 

En la actualidad existen diversos modelos 
físico-matemáticos basados en ecuaciones 
diferenciales o algebraicas, que permiten la 
descripción de las cuatro fases del riego por 
melgas. Entre los modelos disponibles se 
puede mencionar el RIGRAV (Rendón et al., 
1997), que utiliza una combinación del modelo 
hidrológico para describir el flujo superficial, 
y la ecuación de Green y Ampt (1911), para 
modelar el flujo del agua en el suelo; otro 
modelo disponible es el BRDFLW (Strelkoff, 
1985), que utiliza la forma completa de las 
ecuaciones de Saint-Venant, para representar el 
flujo superficial, y una ley de infiltración tipo 
Kostiakov, para representar el flujo del agua 
en el suelo, es decir, una forma que no tiene 
base físico-matemática para la descripción 
del fenómeno. El modelo desarrollado por los 
autores del presente trabajo tiene utilidad con 
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fines de diseño del riego por melgas, con la 
ventaja de incorporar dos ecuaciones con base 
físico-matemática para describir el proceso: las 
ecuaciones de Saint-Venant, para describir el 
flujo sobre la superficie del suelo, y la ecuación 
de Richards, para modelar el flujo del agua en 
el suelo. 


Eficiencias en el riego por melgas 


En riego es indispensable distinguir al menos 
tres eficiencias relacionadas con el proceso, 
la eficiencia de aplicación, la eficiencia de 
requerimiento de riego y la eficiencia de 
La eficiencia de 
aplicación (Ea) se define como Ea = Vr/Vp, 
donde Vr es el volumen requerido para satis- 


uniformidad del riego. 


facer las necesidades de agua en la zona de 
raíces del cultivo [L*] y Vp es el volumen de 
proyecto [L*]. Vr se obtiene con la expresión: 
Vr = LnAr, donde Ln es la lámina de riego neta 
[L], definida conforme a los requerimientos de 
riego del cultivo, y Ar es la superficie de riego 
considerada [L?]. El volumen de proyecto se 
calcula como Vp = QpTr = LbAr, donde Qp es 
el gasto de proyecto [L*T*], Tr es el tiempo de 
riego [T] y Lb es la lámina de riego bruta. Si en la 
expresión (18) el numerador y el denominador 
se dividen entre Ar, se obtiene Ea = Ln/Lb. 

La eficiencia de requerimiento de riego 
(Er) se define como Er = Vd/Vr, donde Vd es el 
volumen disponible [L*]. Esta eficiencia indica 
la manera en que se están satisfaciendo las 
necesidades de agua del cultivo. 

La situación ideal respecto a la uniformidad 
se presenta cuando todas las plantas reciben 
la misma cantidad de agua, situación que 
equivale a aplicar una lámina de riego uni- 
forme en toda la longitud de la melga. Para 
evaluar la uniformidad en la distribución de la 
lámina infiltrada puede utilizarse el coeficiente 
de uniformidad de Christiansen (CUC): 


(18) 


donde 1, es la lámina infiltrada en el punto 1 
[L]; Í es la lámina infiltrada media [L], y nes el 
número de puntos considerados para efectuar 
el cálculo. 


Cálculo del gasto óptimo 


Empleando un modelo numérico desarrollado 
para el acoplamiento de las ecuaciones de 
Saint-Venant y Richards, es posible determinar 
el gasto para el cual se obtiene la mayor efi- 
ciencia de uniformidad, manteniendo valores 
lo más elevados posibles de las eficiencias 
de aplicación y de requerimiento de riego; 
es decir, para determinar el gasto óptimo de 
riego, teniendo en cuenta diferentes láminas 
de riego, longitudes de melga y tipos de suelo. 
En una primera etapa se realiza la estimación 
del gasto óptimo para el suelo franco de 
Montecillo reportado en la literatura (Fuentes 
1992); en este caso se tienen parámetros para la 
ley de resistencia de Fuentes et al. (2004), d = 1 
(régimen de flujo laminar), k =1/54; parámetro 
en la ecuación de cantidad de movimiento, f 
= 2; parámetros para la característica de van 
Genuchten (1980), 0, = 0.4865, 8, = 0.0, m = 
0.126, n = 2.288, y, =-32.75; parámetros para la 
conductividad de Brooks y Corey (1964), K, = 
1.84 cm/h, n = 11.0; valor inicial del contenido 
volumétrico de agua 0, = 0.2749; pendiente 
topográfica J, = 0.002 cm/cm. 

Se ha utilizado un paso de tiempo constante 
At=1.0s para el acoplamiento de las ecuaciones 
de Saint-Venant y Richards. En el caso de esta 
última, se utilizaron pasos de espacio mínimo 
y máximo: Az_. =0.02 cm y Az,. =1.0 cm. La 
discretización de la ecuación de Richards se ha 


mín máx 


efectuado haciendo uso de valores de pasos de 
tiempo y espacio, que generan una malla más 
fina que la utilizada por Saucedo et al. (2002), 
con lo cual se espera mantener una cota de 
error relativo máximo menor que el 0.25%, 
reportado por dichos autores. La discretiza- 
ción utilizada para la solución de la forma 
completa de las ecuaciones de Saint-Venant 
guarda semejanza con las reportadas en la 
literatura: Katopodes y Strelkoff (1977), At. 
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= 5 s; Akanbi y Katopodes (1988), Af_. =1 s; 
Playán et al. (1994), At__ =2.12s. 

En la figura 2 se muestran perfiles de 
flujo típicos para la fase de avance obtenidos 
mediante el modelo desarrollado por los au- 
tores de este trabajo; en la figura 3 pueden 
apreciarse curvas de avance y recesión típicas; 
en la figura 4 se muestra un detalle de la 
evolución de tirantes al final de la melga 
durante la fase de almacenamiento. 


Relación entre la longitud de la melga y el 
gasto óptimo 


La eficiencia de uniformidad, medida a través 
del coeficiente de uniformidad de Christiansen 
(CUC), puede obtenerse para diferentes com- 
binaciones de longitud y gasto de aporte en la 
cabecera de la melga. En la figura 5 se muestra 
un ejemplo de los valores que se obtienen para 
cuatro longitudes de melga; en el caso del suelo 
franco de Montecillo, puede apreciarse que la 


eficiencia de uniformidad varía de manera 
sensible con el gasto de riego. 

Para cada longitud de melga es posible 
determinar el valor del gasto de aporte que 
produce un máximo en el coeficiente de 
uniformidad, manteniendo valores lo más 
elevados posibles de las eficiencias de apli- 
cación y de requerimiento de riego; el gasto 
así determinado se denomina gasto óptimo. 
Cuando se hace variar el gasto de aporte, las 
eficiencias de aplicación y de requerimiento de 
riego no varían de forma significativa, no así 
la eficiencia de uniformidad, la cual cambia 
sensiblemente con el gasto de riego; es decir, las 
eficiencias de aplicación y de requerimiento de 
riego puede considerarse que no son variables 
de decisión en la definición del gasto óptimo y, 
por lo tanto, la eficiencia de uniformidad es la 
que permite definir el gasto óptimo de riego. 

La cantidad de agua que es proporcionada 
en la cabecera de la melga puede calcularse 
como el producto del gasto de aporte (q,) por el 
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Figura 2. Perfiles de flujo obtenidos durante la fase de avance, correspondientes a tiempos de 1, 5, 10 y 20 minutos. 
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Figura 3. Curvas de avance y recesión. 
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Figura 4. Detalle de la evolución de tirantes al final de la melga durante la fase de almacenamiento. 
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Figura 5. Variación de la eficiencia de uniformidad para diferentes longitudes de melga y gastos de aporte 


en el suelo franco de Montecillo. Lámina de riego de 10 cm. 


tiempo de riego (1,), mientras que la cantidad 
de agua que ha sido aplicada en la melga es 
proporcionada por el producto de la longitud 
de la melga (L) y la lámina de riego (1); el 
balance de la masa proporciona: 


lo T, =Lf, (19) 


entonces: 
(20) 


si se especifican el tipo de suelo, la pendiente 
topográfica de la melga y el factor de fricción 
de la ley de resistencia hidráulica, es posible 
encontrar una relación lineal entre el gasto de 
aporte y la longitud de la melga, es decir: 

q, =CL (21) 
la comparación de las ecuaciones (20) y (21) 
implica que: 


(22) 


como c es un valor constante, se deduce que 
para la aplicación de una lámina de riego 
específica existe un tiempo de riego único para 
obtener un valor máximo del coeficiente de 
uniformidad. 

Al correlacionar diversos valores de gastos 
de aporte óptimos y longitudes de melga, es 
posible establecer que existe una proporción 
básicamente lineal entre ambas variables para 
un suelo homogéneo. En la figura 6 se muestra 
la relación que guardan el gasto óptimo y la 
longitud de la melga para el suelo franco de 
Montecillo; las láminas aplicadas son de 8, 10 
y 12 cm, debe apreciarse que existe monotonía, 
en el sentido de que la pendiente de la relación 
entre la longitud de la melga y el gasto óptimo 
disminuye conforme se incrementa la lámina 
de riego. 


Cuadro de diseño del riego por melgas 
Es posible obtener relaciones entre el gasto 


óptimo y la longitud de melga para diversos 
tipos de suelo, diferenciados según su ca- 
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Figura 6. Relación entre la longitud de melga y el gasto de aporte óptimo para el suelo franco de Montecillo 


para tres láminas de aplicación: 8, 10 y 12 cm; Ks en cm/s. 


racterización hidrodinámica, procediendo 
como a continuación se indica. El contenido 
volumétrico de agua residual ha sido asumido 
igual a cero de acuerdo con Fuentes et al. (1992). 
El contenido volumétrico de agua a saturación 
ha sido asimilado a la porosidad total del 
suelo ((), determinada con base en los valores 
proporcionados por Rawls y  Brakensiek 
(1981), en relación con la textura del suelo. 
Para determinar el valor de los parámetros m 
y n, se reconstruyó una curva granulométrica 
para cada suelo, con base en los porcentajes de 
arena, limo y arcilla presentes en el triángulo de 
texturas, y se siguió el procedimiento sugerido 
por Fuentes (1992). La conductividad hidráu- 
lica a saturación, así como la escala de presiones 
Y, estimada en relación con el parámetro de 
succión en el frente de humedecimiento (h y) de 
la ecuación de Green y Ampt, fueron estimados 
para cada suelo mediante los gráficos de Rawls 
y Brakensiek (1981). 

El contenido volumétrico de agua necesario 
para asignar la condición inicial para la ecua- 
ción de Richards se determinó teniendo en 


cuenta la humedad aprovechable de cada tipo 
de suelo, suponiendo que se ha consumido el 
50% de dicha humedad antes de la aplicación 
del riego. La humedad aprovechable se deter- 
minó restando los contenidos volumétricos 
de agua correspondiente a la capacidad de 
campo (CC) y al punto de marchitamiento 
permanente (PMP); los valores de CC Y PMP 
para cada tipo de suelo se calcularon conforme 
alos reportados por Rawls y Brakensiek (1981), 
asociados con el triángulo de texturas. 

En el cuadro 1 se muestran los valores de 
los parámetros involucrados en las carac- 
terísticas hidrodinámicas para diez tipos de 
suelo. El conocimiento de los valores de los 
parámetros involucrados en las características 
hidrodinámicas para cada tipo de suelo y del 
valor del contenido volumétrico de agua inicial 
hace posible calcular las relaciones entre la 
longitud de la melga y el gasto de riego óptimo 
para cada suelo, dada una lámina de riego; 
el valor de la constante que relaciona ambas 
variables se expresa en términos de gasto de 
aporte por unidad de área, es decir, por unidad 
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de ancho y por unidad de longitud de melga; el 
resultado se muestra en el cuadro 2. 

Rendón et al. (1997) presentan un cuadro 
de diseño similar al cuadro 2; sin embargo, el 
cuadro reportado en la literatura tiene algunas 
inconsistencias de monotonía entre la relación 
que guardan las variables gasto óptimo, tiempo 
de riego y lámina aplicada; si bien son un 
porcentaje pequeño, los resultados que rompen 
la monotonía pueden reflejar que el modelo 
utilizado para la elaboración del cuadro de 
diseño del manual presenta dificultades para 
reproducir el fenómeno para tiempos de riego 
grandes. El acoplamiento de las ecuaciones 


de Saint-Venant y Richards permite obtener 
resultados que guardan la monotonía en las 
variables del diseño, según puede apreciarse en 
el cuadro 2. En el mismo sentido, se observan 
diferencias en cuanto a los gastos de aporte y 
los tiempos de riego obtenidos en este trabajo 
y los reportados en el manual citado; dichas 
diferencias son más significativas cuanto más 
arcilloso es el suelo, situación que se tenía 
prevista debido a que la ecuación de Green y 
Ampt utilizada en el modelo empleado para 
desarrollar el cuadro de diseño del manual 
representa con mayor suficiencia el flujo de 
agua en suelos arenosos. 


Cuadro 1. Características hidrodinámicas de suelos para el diseño del riego por melgas. 


Textura del suelo 8, (cm*/cm?) 9, (cm?/cm*) Ks (cm*/h) n m y, (cm) 
Arcilla 0.525 0 0.010 61.17 0.0229 132.50 
Arcilla limosa 0.500 0 0.015 SO 0.0440 94.70 
Franco-arcillo-limoso 0.500 0 0.070 15.33 0.0905 57.80 
Franco-arcilloso 0.475 0 0.150 19.30 0.0714 34.15 
Arcilla arenosa 0.425 0 0.200 41.47 0.0327 23.70 
Limo 0.500 0 0.500 5.60 0.2477 30.70 
Franco 0.475 0 0.700 13.93 0.0989 19.20 
Franco limoso 0.525 0 0.600 12.01 0.1165 29:39 
Franco-arcillo-arenoso 0.425 0 1.500 18.44 0.0736 33.99 
Franco arenoso 0.450 0 5.000 13.62 0.1004 52. 


Cuadro 2. Cuadro de diseño del riego por melgas: gasto en 1/s/m? para la aplicación Óptima de la lámina de riego. 


Factor de fricción k= 1/54, pendiente topográfica J, = 0.002. 


Lr=8 cm Lr=10 cm Lr=12 cm 

Textura del suelo a, Tr a, Tr a, Tr 

(1/s/m?) (h) (1/s/m?) (h) (1/s/m?) (h) 
Arcilla 0.00012 224.1 0.00010 338.2 0.00009 445.0 
Arcilla limosa 0.00014 201.6 0.00012 270.5 0.00011 Sono 
Franco-arcillo-limoso 0.00060 44.1 0.00050 66.6 0.00046 82.9 
Franco-arcilloso 0.00088 31.4 0.00078 44.0 0.00072 57.8 
Arcilla arenosa 0.00090 28.7 0.00080 42.4 0.00077 52.0 
Limo 0.00399 6.9 0.00333 10.0 0.00296 Sí 
Franco 0.00411 6.4 0.00354 9.6 0.00326 125 
Franco-limoso 0.00446 6.2 0.00388 8.8 0.00349 11.6 
Franco-arcillo-arenoso 0.00490 5.8 0.00476 7.4 0.00464 9.0 
Franco-arenoso 0.02476 1.2 0.02223 1.6 0.02073 2.0 
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Conclusiones 


Se verificó la relación que existe entre la longi- 
tud de una melga y el gasto óptimo de riego que 
debe aplicarse para obtener un valor máximo 
en el coeficiente de uniformidad, así como 
valores elevados de la eficiencia de aplicación 
y de la eficiencia de requerimiento de riego de 
los cultivos. La forma lineal de la proporción ha 
sido reportada en la literatura, como producto 
de la aplicación de un modelo hidrológico 
para el flujo a superficie libre, y la ecuación 
de Green y Ampt, para describir el flujo del 
agua en el suelo; en este trabajo se verifica la 
relación mediante un modelo hidrodinámico 
completo, que hace uso de las ecuaciones de 
Saint-Venant para la descripción del flujo del 
agua sobre el suelo acopladas internamente 
con la ecuación de Richards para la modelación 
del flujo del agua en el suelo. Los coeficientes 
de proporcionalidad calculados difieren de 
forma significativa, según el modelo utilizado, 
debido básicamente a la diferencia entre la 
mayor complejidad del modelo empleado 
en este trabajo y la del modelo aplicado en 
la literatura. De acuerdo con los resultados 
obtenidos mediante la aplicación del modelo 
para diez tipos de suelo diferentes, se incluyen 
cuadros de diseño operativo del riego por 
melgas. 
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Abstract 


SAUCEDO, H., ZAVALA, M., FUENTES, C. é: CASTANEDO, V. Optimal flow model for 
plot irrigation. Water Technology and Sciences (in Spanish). Vol. IV, No. 3, July-August, 
2013, pp. 135-147. 


The primarily linear relationship was verified between the length of a plot and optimal 
irrigation flow, which is the flow needed to obtain the maximum value of the uniformity 
coefficient, resulting in high values for application efficiency and the amount of irrigation 
required. This verification was performed by applying a complete hydrodynamic model using 
numerical methods, based on Saint-Venant equations to describe water surface flow and 
internally coupled with the Richards equation in order to model it. Operating design tables 
for plot irrigation were developed based on the results obtained from applying the model to 10 
different types of soil. 


Keywords: plot irrigation, Richards equation, Saint-Venant equations, optimal irrigation 
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